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Abstract — The flow field and mixing characters of high 
efficiency vortices in one type of Radially Stratified Flame 
Core Burner RSFC where analysed through Cold 
computing fluid dynamics model CFD in accordance with 
Mathematical model for flame visualization. The model 
gives good results and clears the important aerodynamics 
features of a typical low Nitrogen Oxides NOx of internal 
staging schematic. The model can be used to evaluate this 
design performance or to get the necessary information for 
development its efficiency in accordance with the global 
demand to clean energy with high efficiency and low NO x 
emissions. The mixing of this burner dominated by coherent 
structures, so it is outperform the molecular diffusivity and 
the length which fluids will mix is independent of the fluid 
velocity and can determined only from Taylor macroscale 
which related by the physical size of the fluid domain or in 
better understanding with the internal recirculating zone 
dimensions because it is form a region for axial mixing for 
species and smear out their residence time distribution. 
Keywords — coherent structure, gaseous mixing 
characters, internal staging schematic flame, cold flow 
flame model, Taylor macroscale, internal recirculation 
zone, penetration theory. 

I. INTRODUCTION 

In the present work, it appears that a statistical analysis of 
the mixture quality has been sufficient to obtain information 
which can be used in a relatively fuel/ air flow adjustment 
scheme to minimize the NOx emissions. Swirl burners in 
industrial furnaces utilize powerful vortices to increase the 
speed of collision between axial and tangential flows, thus 
speeding up the time for mixing fuel and air and also 
extending the residence time. Hence, as suggested by Coats 
(Coats et al. 1996) the mixing of systems dominated by 
coherent structures cannot be modelled accurately with 
gradient based techniques. As in most technical 
applications, the stirring effect of large structures will 


namely outperform the molecular diffusivity. The problem 
can become even more complex in the case of counter 
gradient diffusion processes which occur during combustion 
(Veynante et al. 1996). 


II. METHODOLOGY 



Fig.l: typical low NOx RSFC flow field with internal 
staging schematic [3] 

The mathematical model for fluid dynamics field shows one 
type of Radially Stratified Flame Core Burner RSFC which 
clears the important aerodynamics features of a typical low 
Nitrogen Oxides NOx of internal staging schematic see 
fig(l).The burner used in this study consists of two 
concentric tube with a tangential and radial swirler placed in 
annular arrangement .So an initial fuel -rich flame core is 
created by the mixing process between the central fuel jets 
with the primary air which is a portion of air combustion 
equals to 30% tends to move around the air deflector as 
primary air , and 70% of the air combustion enters through 
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a radially displaced annulus as secondary air. The fuel jet 
penetrates through the internal recirculation zone as seen in 
fig. (2) which clear that there are two mixing regions for the 
jet in crossflow, the first is the jet momentum dominated 
near field (MDNF), and the second is the main flow 
momentum dominated far field (MDFF ) as it was 
mentioned by lee and Chu et al 2003 who also introduced 
the length scale l mv to separate the two zones: 

^mv ~ JpJ Mvo / Pa u a 

Where M vo = uf sin(0 ; ) n df /4 is the vertical 

momentum flux of the jet at y = 0, u a the bulk velocity of 
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the main flow, pj and p a are the densities in the case of two 
fluids. For the momentum dominated near field the jet 
penetration y verifies y ll mv < 1, and for the momentum- 
dominated far field, y ll mv > 1. In a similar way, they define 
an integral length scale for a jet in co-flow flow, 

— JpjMjo/PaMa 

Where M J0 = Uj(uj — u a ) sin(0y) n df /4 and they define 
for the momentum dominated near field (MDNF) the jet 
penetration y verifies x/l ^ < 10 is valid. 



Fig. 2: mixing regions for the jet A [4] and structures which are generated by a jet in crossflow B [4], [5]. 


Shortage of detailed and accurate experimental data on fuel-air mixing in furnaces is due to the difficulty and complexity of 
measurements in flames. 

The swirl of fluids in the low NO x burner forces more mixing process. So its effects is thank to create turbulent operating 
conditions for mixing fuel /air by inducing high efficiency vortices in addition to benefit from the blades which twist the fluid 
away or toward in helical movement depends on its position when it flow axially toward the burner 

For sufficient turbulence and mixing, the oxidizer try to enter into the fuel rich region by laminar or turbulent diffusion at the 
boundary of the fuel and oxidizer jet or it take another way in the field of internal recirculation zone which recycle flue gas back 
into the flame (see fig.(3)) 
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Recirculation caused by wake behind a bluff body 



Axial confined jet and secondary recirculation 
Fig. 3: In internal FGR, the fuel gas is recycled into the 
flame zone due to burner aerodynamics to reduce NO x [6]. 

In all combustion systems, according to Zeldovich 
mechanism it is usual to have little quantities of thermal 
NOx emission at flame temperature which is below 1450 
°C, and that is because the oxygen atoms will not be 
available from their molecules to start thermal NOx 
formation so the combustor parameters like the residence 
time, flame temperature, pressure and heat losses which are 
relevant in the generation of NOx must be adjusted for the 
global demand to clean energy with high efficiency and low 
emissions. 

A minimization of NOx emissions in combustors is 
achieved when the fuel/air mixing quality is optimized, as 
homogeneous as possible but this is not true for mixtures at 
stoichiometric conditions [7] (see fig.(4 )) 



equivalence ratio 0 



0.4 0.8 1 1.2 0 


Fig. 4: impact of fuel/air mixing quality on the NOx 
emissions for different fuel/air equivalence ratios cp, from 
Joos B [8] and Dominant NOx formation mechanisms 
shows Five chemical mechanisms describe how NOx are 
produced in combustors A [8]: *thermal NOx (or 
Zeldovich) *Fenimore Mechanism (also called prompt NOx 
) *Fuel NO by N 2 Oroutes *NNF[ Radical 


III. MATHEMATICAL FLAME 
VISUALIZATION MODEL 

In this combustor the presence of the internal recirculation 
zone plays important role in flame stabilization and 
introducing rotation in the stream by the vanes which are 
mounted on a central gaseous fuel tube by the air deflector 
geometry has a big role in mixing process in the middle of 
burner, this effect is done by controlling the aerodynamics 
of the internal recirculating zone to mix fresh fuel with the 
primary combustion air. But the vanes which occupy the 
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space in the annular around it by a vane type - swirler are 
thank for playing a mazing effects on recycle the secondary 
combustion air into reaction zone .( see fig. (5))Which 
shows the geometry of the burner. A swirl burner have 
convergent-divergent air exit shaped section (also called 
quaral or burner throat section) , the air from the swirler is 
contracted first in the air exit and then it spreads into the 
furnace by using the expansion cone . 



Fig. 5: Shows a schematic diagram of the air and fuel flow 
velocities of a swirled jet burner (a) [9], burner geometry 
(b), and (c) a proto type of swirl burner, its components and 
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Confusor, which for similarity conditions the all were done 
in line with fact. 

The volumetric stoichiometric fuel/air mixture ratio is 
computed as [9] 



A nblu c 

Where u F the fuel velocity from the nozzle with diameter 
d F and the air mixing with fuel stream is characteristic by 
mixing velocity u c within a cylindrical volume its length 
equal the flame length l and its wide equal b which is the 
boundary where the average axial velocity is zero. 

The characteristic mixing velocity u c is the sum of two 
items, one is due to recirculation and another item is due to 
shearing between fuel and air stream. So it can be computed 
as [9]: 

u c = n. b. 1. u RZ + constant. d F . 1. n. ( u F — u A ) 
Where u RZ the velocity of characteristics recirculation zone 
and is computed from the relation [9]: 

rb 

u RZ = I —u z .2n.r.dr/nb 2 

Jo 

Where u z the mean axial velocity and b is located where 
the flow is perpendicular to the jet axis. So by replace the 
u c by its two items in the relation which computed the 
mixture ratio, it can be obtained the following one [9] 



d A f u rz ■ ft _|_ (Ma ~ u f) d F ^ \ 

\u A .d A u A d A 2 ) 

Where subscript F, A denotes values for fuel and air stream, 
m, d, u, 1 are mass flow, diameter, velocity respectively and 
1 flame length while c t , c 2 are experimental constants. 

So from the last equation it was seen that burners’ geometry, 
viscous shearing, swirl and kind of fuel all of them affect 
the flame features and their effects can be seen by applying 
computing fluid dynamics model to notice the flame 
structure s’effects of on NO x emissions. 

So the following are the technical data of burner: 

- Burner axis: vertical 

- Regulating range of capacity of burner 20%-60% 

- Excess of air nair A =1.15 

- Required type of the burner: gas type with natural draft of 
air fitted with a gas injector stabilization torch, Hydrogen 
content in the heating gas must not exceed 60% volume, 
and the fuel doesn’t contain any nitrogen compounds 

- Chemical composition %-vol for gaseous fuel: H2-28, Ci- 
34, C 2 -17, C 3 -17, 03, C 5 - 1, H 2 0-1.5, H 2 Smax.5mg/Nm 3 

* Net calorific value: llOOOkcal/kg 
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* Density kg/Nm 3 0.924 

* Flammability limits% 3.4-16.5 

- Temperature in the spot of the burner inlet 30°C 

- Heating gas consumption of I burner NrnVh; 

For min. output 23.2, for nominal output 58.0 and for max. 
Output 69.6 

IV. COLD COMPUTING FLUID DYNAMICS 
MODEL CFD 

- The interface F is between two fluids (air /fuel) but it is 

movable and unknown. So, it must determine by the 
solution of CFD model (see fig. (6)). 



Fig. 6: Schematic of a two-fluid domain [10] 


- As in the case of any problem in flow filed this particular 
type of flow described by the mathematical expressions of 
continuity and momentum in Cartesian coordinates 
The continuity equation: 


— + V. ( pv ) = 0 

OT 

Where p is mass density for each fluid (fuel/air), r is the 
time and ~v = ( u,v,w ) is the three velocity field, and by 
writing the local component yields by using the mass 
fraction w for each fluid [11] 

= -v./ 

Where J is the Fickian diffusive mass flux of fluid 1 in 
corresponding fluid: 

/ = —pD. Vw 



As the mixing process is assumed isothermal, the density is 
a function of pressure p and composition, so using the 
derivatives of density in continuity equation yields [11]: 
dp tdp \ dp (dw ^ \ 

a? (si + ”• Vp ) + (aF + ”• Vw ) + ”■ v ' v = 0 


Neglecting the spatial variation of p on the density and 

rearranging the last equation by using the equation of spatial 

variation of w for the fluid 1 yields [11] 

_ -1 dp dp 1 dp 

V.v = —. — . — + —V./ 

p dp dt p 2 dw 

This equation indicates that compressibility and non¬ 
ideality (density variations with compositions) can both be 
source of divergence of the velocity field [11] 

For Computing cold fluid dynamics model some 
assumptions are used as assume that the densities of the 

fluids are uniform at each side of the interface P so that the 


Atwood number At is very small: 


At = (Pi - p 2 )/(Pi + p 2 ) « 1 
So the Bossinsesq approximation can be made in this case 
[10], and it retains density variations in the gravitational 
term (giving buoyancy forces) but disregards them in the 
inertial term; i.e. in the vertical momentum equation: 


du t d(U]U^) 
dr dXj 


dp d T u 


Where: T is the actual temperature applicable locally. : 
is the undisturbed constant temperature and by using the air 
as working incompressible perfect gas it is possible to write 
B t : is the coefficient of the thermal expansion of fluid and 
by using the equation of stat for the perfect gas 

V 

P ~ RT 

Where R is universal gas constant, p is the local static 
pressure, and p is the local density so ( B T = 1 /T), T t j is 
the Reynolds stress tensor divided by the density and it is 
computed by using the Boussinesq model equation 


T ij = -u i u j =v T — + 


du t i dUj 


dxj dx t 3 dx { 


2 du k \ 2 

k 8u)--KS u 


Where K is the turbulent kinetic energy, 8 t j is Kronecker 

Delta Function, Then the letters i, j, k ...can be used as 
variables, running from 1 to 3 so Instead of using x, y, and z 
to label the components of a vector, we use 1, 2, 3 and v x is 
the turbulent eddy viscosity can be calculated by K-Epsilon 
model so 

2 

K 

v T -C^ — 

Where C^ is constant equals to 0.09 
These models are not valid near wall to model wall effect so 
wall function have to be employed. 

4-1 wall y+ strategy for dealing with wall-bounded - 
turbulent flows: 
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It was observed that viscous sublayer turbulence model 
were in better accord with results while the model results 
obtained from the wall function bridge the gap with the 
walls were not suitable choice for the internal recirculation 
field simulation and this is related by the effect of positive 
and adverse pressure gradient in the direction of flow [12]. 

It is also advisable to avoid having the wall-adjacent mesh 
in the buffer region since neither wall functions nor near 
wall modelling approach accounts for it accurately, so for 
the existence of deferent boundary layers in burner 
geometry an enhanced wall function is used. 

4-2 Boundary conditions: Boundary conditions are used 
according to the need of the model. The inlet velocities are 
used similar to the max fuel gas flow conditions in order to 
have a comparison. General correlations are used to 
estimate the turbulence boundary conditions which are 
specified by estimating the turbulence intensity and length 
scale (I = 0.16Re“ 1/8 ) and (1 = 0.07L). 

V. RESULTS AND DISCUSSION 

Velocity vector plots can be seen below in Fig. 8 (a), (b) 
and (c). These plots give an idea of flow separation at 
region which surrounds the burner tip in the middle of 
burner. The recirculation at the wake region of air deflector 
is also obvious from them, the iso surface plots of velocity 
give the most important huge quantities information about 
the movement of fluid were its axial velocity is zero in 
another expression where it is possible to try turn down to 
the reaction zones (flame) (see Fig. 8 (a)) 


The boundary of the recirculation eddy is determined by the 
radial points at which the forward mass flow equals the 
reverse mass flow at the axial station so the region which is 
enclosed where the axial velocity reaches zero magnitude in 
each axial plane [13] (see fig 7) .So The internal staging 
schematic adjusts the flame flow field to control the NOx 
emission that is because the high swirl provides an intense 
internal recirculating zone which recycles flue gas back into 
the flame and thus reducing the oxygen concentration as 
well as lowering flame temperature to prevent the NOx 
emission. 

The simulation was done in steady frame by using NX 
Siemens code which depends near wall modeling 
assumptions for each turbulent model because some 
turbulent cases needs fine mesh close to the walls to capture 
the flow characteristics field as in the case of presence flow 
separation and/or not fully- developed flows. So enhanced 
wall function where applied on the walls of burner 
components especially for nozzles, burner tip and air 
deflector walls where to have accurate flow field model 
their wall-adjacent mesh obeyed this condition y+ < 5 or y+ 
>20 (see Fig. 8 (b)) 

Images in fig. 8 (b) show that as the tangential momentum is 
added to flow, radial pressure gradients increase. Essentially 
the fluid wants to move out from the center of rotation, 
these radial pressure gradients generate a low axial pressure 
zone which then draws material back up into a CRZ. 



Fig. 7: Sketch of the recirculation zone in a confined 
combusting flow. 
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Load Case 1, Static Step 1 

Velocity - Element-Nodal, Averaged, M, 

Min : 0.01, Max : 240.10, Units ■ m/sec 
- 240.10 


Load Case 1, Static Step 1 

Velocity - Element-Nodal, Averaged, Magnitude 

Min : 0.01, Max : 240.10, Units - m/sec 


Maximum 
Nude 1000940 
210.D97 nV99C 


Minimum 
Mode 355672 
|o DU m/sec 


: Solution 


as*embJy1_fem2_sim1 sim last one : Solution 5 Resu assembly1_fem2_sim1.sim last one : Solution 5 Resu 


assembly 1_fem2_slm1 .slm last one : Solution 5 Result 
Load Case 1, Static Step 1 
Velocity - Element-Nodal 
Min: 0.01, Max: 240 10, 


216.&9 
192.08 
168.07 
144.06 
120.05 
96.05 

■ 7240* - ^ 

48.03 
- 24.02 
" - 0.01 

Units ® m/sec 
as sembly I Jem2_sim 1. s im 
Load Case 1, Static Step 
Velocity - 

Min: 0.01, Max: 240. 


n 

Load Case 1, 
Velocity - Ele 


c ' 


assembiyl Jem?_;im1 si 
LoadCa^s 1. Static Step 
Velocity - 
Min: 0.01, Max 


Fig. 8: a plots of isothermal surface of velocity 



assernbly1_fem2_sirn1 .sim last one : Solution 5 Resi 
Load Case 1,. Static Step 4 
Y Plus - Element-Nodal, AverSgeiJJScalar 

Min : 0.114 Max : 2.13E+004, Units = 

■ 2.13E+0D4 M ^ 

4 22E+003 fl 

i:88E+003 ^ l» 


! 1, Static Step 1) 
Element-Nodal.llj 
EVW : 239 ,54$l 




assembly 1_fem2_sim1.sim last one : Solution 5 Result 
Load Case 1, Static Step 1 
Velocity - Element-NodaLJfc^i&ged, Magnitude 
Min : 0.01, Max : 

JSfSK*?*?/ 


Solution 5 Result 




assembly1_fem2_sim1 .sim last one : Solution 5 Resu 

Load Case 1, Static Step 1 

Velocity - Element-Nodal, Averaged, Magnitude 

Min : 0.012, Max : 240.10, Units = m/sec 

1 240.10 
124.34 
64.40 
33.35 
17.27 
8.94 
4.63 
2.40 
1.24 
0.643 
0.333 
0.173 
0.089 
, 0.046 
0.024 
1 0.012 
Units = m/sec 


assembly 1_fem2_sim1 .sim last one : Solution 5 Result 
Load Case 1, Static 
Velocity- 
Min: 0.01,1 


Fig. 8: (b) shows the recycle of fluid due to burner aerodynamics to reduce NO x in spite of the mesh quality and y + . 


It is clear that fuel velocity exceeds laminar flame, and an ignition ledge is used as air deflector and even if the air flows at very 
high speed at maximum fuel rate , the air speed very close to the ledge will be small enough , so the flame can establish very 
close to the ledge ignition and be quit stable even over a wide range of firing rate .in other words the flame speeds and flow 
velocity are matched in this recirculating zone and a good mixing of reactants with air occurs due to high turbulence generated by 
high shear stresses ( see fig. 8 (c)) 
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assembly 1_fem2_sim1 .sim last one : Solution 5 Result 
Load Case 1, Static Step 1 

Turbulence Dissipation - Element-Nodal, Averaged, Scalar 
Min: 0.89, Max 

I 2.55E+008, 

3.64E+00, 

- 5.19E+0J 
1 7.41 E+j 
1.061 

- 1.51E+1 
* 2.15E+003 
! 3.08E+002 
! 43.89 

I 6.26 

■ 0.89 

Units = m A 2/sec 


assembly 1_fem2_sim1. sim last one : Solution 5 Result 
Load Case 1, Static Step 1 

Turbulence Energy - Element-Nodal, Averaged, Scalar 


assembly 1_fem2_sim1 .sim last one : Solution 5 Result 
Load Case 1, Static Step 1 
Velocity - Element-Nodal, Ave 



1 0.012 

Units = m/sec 


Fig. 8: (c) shows that flame speeds and flow velocity are matched in the recirculating zone (note plots on two paths) 


The flow fields of jets in crossflow are quite complex and fully three-dimensional (full turbulent ), as shown in Fig.(9) They 
involve four main features [5], [14]: "(1) the horseshoe vortex, due to the adverse pressure gradient in front of the injection hole, 
(2) the wake vortices, (3) the shear layer ring vortices, which is induced by a Kelvin-Helmholtz instability in the shear layer of 
the jet, and (4) the counter-rotating vortex pair (CVP), generated by the deflection of the jet, and which seems to be induced by 
the shear layer vortices. The latter structure dominates the far field mixing, even if it is generated in the near field of the jet" [7]. 
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I assembly 1_fem2_sim1.sim last one : Solution 5 Result 
I Load Case 1, Static Step 1 
| Vorticity - Element-Nodal, Averaged, Magnitude 

Unitless ^ 


assembly 1 _fem2_sim 1 .shjvdaSf one : Solution 5 Result 
Load Case 1, Static Step 1 
Vorticity - Element-Nodal, Averaged, Magnitude 
Min : 0.045, Max : 6.01 E+005, Units = Unitless 
Ctgftrt BflAX^york Rectangular 
1 1 . 17E+005 
■ 2.26E+004 
4.38E+003^- 
3.49E'+002 
164.63 

31.92 ~V- 1 

6.19 . \ 


assembly 1_fem2_sim1 .sim last one : Solution 5 Result 
Load Case 1, Static Step 1 

Y Plus - Element-Nodal, Averaged, Scalar - 

Min : 0.114, Max^2rT3^+b04, Units = Unitless 


Min : 0.045, Max : 6.01E+005, Units 


- 6.01 E+005 

- 1.17E+005 

- 2.26E+004 

- 4.38E+003 
8.49E+002/3? 

- 164.63 
31.92 Hi 

6.19 

1.20 

- 0.233 

- 0.045 

Units = Unitless 


I Units = Unitless 
assembly 1 _fem2_sirtTtTsini.last or 
L Pad Case 1, Static Step 1 
Velocity - Element-Nodal, Ayfrag 
;Min : 0.012, Max : 240.10, |lk - 
| 240.10 

I 89,48 

" 33.35 
12.43 
4.63 
| 1,73 

' 0.64 j 

f ! 0.240 
| 0.089 
I 0.033 

■ 0.012 /. 
Units = m/sec i 


assembly 1_fem2_sim1.sim last one : Solution 5 Result 
Load Case 1, Static Step 1 
Y Plus - Element-Nodal, Unaveraged, 

Min: 0.114, Max : 2.13E+0( 


lolution 5 Result 


assembly 1_fenY 
Load Cas^djj 

Veloo«BH 


Solution 5 Result 


Unitless 


2.13E+004 
6.32E+fltfif 
1.88E*Q03 

5.58 : E+002 


Units = m/sec 


Fig.9: Shows full turbulent flow in burner and the features of jets flow fields 


The distance from the inlet at which species will be mixed 
l mlx is given by ( l mlx = v. T mix ) where T mix is the mixing 
time and v is the fluid velocity [15] 

( D t oc — and K oc (L r . e) 2/3 ) where K is the turbulent 

energy , L T is the turbulence macro scale, D T is the 
turbulent diffusion coefficient but £ is the turbulent energy 
dissipation rate which is proportional to the multiple of the 
flow rate and pressure drop ,so in result [15] 

( s oc i? 3 ) but ( T mix oc Df 1 ) yields ( £. D T oc L^ /4 . i? 4 ) 

In the result [15] 

L mix oc 1/L 4 / 3 


Thus the length which fluids will mix is independent of the 
fluid velocity and it can determined from Taylor macroscale 
[15] 

The turbulence intensity is larger in regions of hairpin 
vorticity which are generated in high shear regions which 
surrounding the vortices or at the edge of blades (see fig.9 
and fig. 10). The computing Reynolds number is 1.2xl0 5 so 
the flow field is fully turbulent and the turbulent length 
scale will not increase more to the restrictions of physical 
size of the burner so the mixing length in this fully 
established case will remain constant. 
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assembly 1_fem2_sim1.simU^stor»i«Sfili4lign 5 Result 
Load Case 1, Static SteoJ^ 
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Fig. 10: Shows the turbulence intensity of hairpin vorticity which are generated in high shear regions 


VI. CONCLUSION 

Mathematical model studies two aspects of the momentum 
effects on flame which are the forward momentum normally 
associated with the average outlet velocity of the 
combustion products and the lateral momentum caused by 
swirl, So this Mathematical model shows how cold flow 
measurements can be combined with flame visualization to 
model the spatio-temporal response of fuel/air mixing field 
and the ability to use burner aerodynamics as alternative 
design concepts for staging the radially stratified flame core 
burner instead of relying on physically separate zones 

REFERENCES 

[1] Coats, C. M.1996. Coherent structures in combustion. 
Progress in Energy and Combustion Science, 
22(5):427-509. 

[2] Veynante,D. Piana,J. Duclos,J. M. and Martel, C. 
1996. Experimental analysis of flame surface density 
models for premixed turbulent combustion. In 
Symposium International on Combustion, volume 1, 
pages 413-420. Combustion Institute. 

[3] Shihadeh,A.L. Toqan, M.A. Beer, J.M. Lewis P.F, 
Teare, J.D. Jimenez, J.L. and Barta, L. 1994. Low 
NOx emission from aerodynamicall ystaged oil-air 


turbulent diffusion flames. ASME FACT-18, 
Combustion Modeling, Scaling and Air Toxins 

[4] Lee ,J. H. W. and Chu ,V. H. 2003. Turbulent jets and 
plumes - a lagrangian approach. Kluwer Academic 
Publishers. 

[5] Gutmark, E. J. Ibrahim I. M. and Murugappan, S. 
Circular and noncircular subsonic jets in cross flow. 

[6] Fiskum,A. 2008. Calculation of NOx Formation in a 
Swirl Burner. Department of Energy and Process 
Engineering. Norwegian University of Science and 
Technology 

[7] Lapacelle,A.2011. Modeling, control, and 
optimization of fuel/air mixing in a lean premixed 
swirl combustor using fuel staging to reduce pressure 
pulsations and NOx emissions.Berlin university of TU 
,ISBN on line 9783-9783-2382-7 

[8] Johari, H. 2006. Scaling of fully pulsed jets in 
crossflow. AIAA Journal, 44(11):2719-2725. 

[9] Ashworth,D.A.2017.An Analytical Model to Predict 
the Length of Oxygen-Assisted, Swirled, Coal and 
Biomass Flames. Department of Mechanical 
Engineering Brigham Young University 

[10] Lee,H.G. and Kim,J.2011. A comparison study of the 
Boussinesq and the variable density models on 


www.aipublications.com/iireh 


Page | 46 













International journal of Rural Development, Environment and Health Research(IJREH) 
https://dx.doi.Org/10.22161/iireh.2.l.5 

buoyancy-driven flows. J Eng Math (2012) 75:15-27 
DOI 10.1007/sl0665-011-9504-2 

[11] Rongy,L. Haugen,K.B and Firoozabadi, A.2011. 

Mixing from Fickian Diffusion and Natural 
Convection in Binary Non-Equilibrium Fluid Phases. 

Dept, of Chemical Engineering, Mason Laboratory, 

Yale University, New Haven, CT 06511. Wiley Online 
Library, DOI 10.1002/aic. 12685. 

[12] Alabdalah,A.Ibrahim,H.A. Aloosh,E, Hoyas,S.C.2018. 

Wall Y+ Approach for Dealing with Flow into Low 
NOx Swirler Burner in Industrial Furnaces. 

International Journal of Civil, Mechanical and Energy 
Science (IJCMES), Vol-4, Issue-1, Jan-Feb, 2018.[ 

[13] Ashworth,D.A.2017.An Analytical Model to Predict 
the Length of Oxygen-Assisted, Swirled, Coal and 
Biomass Flames. Department of Mechanical 
Engineering Brigham Young University 

[14] Coussement, A. Gicquel, O. Schuller, T and Degrez, 

G.2008. Large eddy simulation of pulsed jet in 
crossflow. AIAA Paper 2010-561, 2010.Physics of 
Fluids, 20:075110, 2008. 

[15] Bakker A., Cathie N., LaRoche R. (1994) Modeling of 
the Flow and Mixing in HEV Static Mixers. 8th 
European Conference on Mixing, September 21-23, 

1994, Cambridge, U.K. IChemE Symposium Series 
No. 136, ISBN 0 85295 329 1, page 533-540 . 


[Vol-2, Issue-1, Jan-Feb, 2018] 
ISSN: 2456-8678 


www.aipublications.com/iireh 


Page | 47 





